Quantum Interference and the Trapped Bose Condensed System 
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In experiments involving Bose condensed atoms trapped in magnetic bottles, plugging the hole in the bottle 
potential with a LASER beam produces a new potential with two minima, and thus a condensate order parameter 
(i.e. wave function) with two maxima. When the trapping potential is removed and the condensate explodes away 
from the trap, the two wave function maxima act as two coherent sources which exhibit amplitude interference. 
A simplified theoretical treatment of this experimental effect is provided by considering momentum distributions. 
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The two fluid Q model has long served as the 
standard picture of superfluidity in liquid He 4 . 
In spite of the great success of the model, Lan- 
dau warned^) that the view was merely quasi- 
classical. It would be desirable to have a quan- 
tum fluid mechanical model in which superposi- 
tion of amplitudes on a macroscopic scale played a 
more central role. For the dilute Bose condensed 
gas[[|, such coherence would mean that the Bose 
condensate order parameter [Q \P(r) would be re- 
quired to maintain many of the properties nor- 
mally attributed to single particle wave functions 
exhibiting amplitude superposition ||. 

Recent experimental probes^, developed for 
observing Bose condensates in atomic traps 
have provided that which is needed to go beyond 
the Landau "quasi-classical" two fluid picture. 
Direct observation of amplitude superposition in 
the Bose condensate order parameter ^(r) has 
been observed^] in mesoscopic samples. 

Consider some central features of recent ex- 
perimental methods. To see what is involved, 
let us first contemplate a "thought experiment". 
Suppose that at times t < CP, there is but one 
atom localized near the origin in the region of 
a trap potential U(r). At time t — 0, the trap 
potential is suddenly removed, and thereafter 
t > + the atomic motion is described by the 
free particle Schrodinger equation ih[dip/dt) = 
-{h 2 /2M)V 2 ip. 

From the viewpoint of quantum mechanics, the 



suddcn[ji(j removal of the trap potential yields 
an out going wave function described using the 
non-relativistic propagator |ll[] G(r ~ s,t). The 
outgoing wave function is determined by 
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From the viewpoint of classical mechanics, the 
particle exits the trap region at a uniform veloc- 
ity along the path r = vt. Both the classical 
(measurement) and quantum (virtual) viewpoints 
are resolved by considering what happens after 
the previously trapped particle freely moves far 
away from the original trapping region. Employ- 
ing |r| >> |s| in Eqs.(l) and (2), one finds that 

■4>{r,t) m (M/{2inht)) 3/2 e lS ( r ' t)/h A(Mr/t), (3) 
where S(r,t) = {M\v\ 2 /2t) and 



A(p) = J e-^ p - s ^ h tP(s,0)d 3 s. 
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is the momentum amplitude of the particle when 
it was trapped by the potential. Eq.(3) implies 
that 



mr,t)\ 2 ^ / \A(p)\ 2 S(r-(pt/M)) 
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Thus, when a particle has escaped away from 
the trapping region (long after the trap poten- 
tial has been removed), the probability density in 
space |-0(r,i)| 2 can be computed by averaging the 
classical equation of motion r = (pt/M) over the 
momentum probability density |^4(p)| 2 , present 
when the particle was in a trapped state. 

To understand recent experiments on quantum 
interference for trapped Bose condensed meso- 
scopic systems, one must replace the above single 
particle amplitudes with their second quantized 
operator counterparts 

d 3 p 
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The many Bosons are held in the trap for times 
t < with an equilibrium momentum distribution 



iv( P ) = (it(p)i(p)>. 
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The atomic interaction potentials have impor- 
tance in assuring that the momentum distribution 
N(p) (for the Bosons within the trap potential) 
actually correspond to thermal equilibrium. The 
experimentalist at time t — suddenly turns off 
the trap potential^^ and optically monitors the 
particle density for times t > 0. The density is 
determined by 



P {r,t) = (ft(r,t)$(r,t)), 
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after the Bosons move away from the original 
trapping potential spatial region. 

To the extent that the interaction potentials 
have a negligible effect on the particles after they 
explode away from the trapping region, the final 
mean density in space p(r, t) is related to the ini- 
tial mean momentum distribution N(p) via the 
many body version of Eq.(5); i.e. 



p(r,*)« / N(p)5(r-(pt/M))- 



d 3 p 
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The ballistic Eq.(9) has been the central theo- 
retical tool employed in the recent measurements 
of trapped Bose condensed systems. The mo- 
mentum distribution of particles within the trap 
is probed by the spatial distribution of particles 
detected away from the trap. 

Let us now review the notion of superposi- 
tion of amplitudes as it becomes apparent in 



the classic "two slit quantum interference exper- 
iment" fl3|| . Quantum particle diffraction ex- 
periments present us with a closely analogous 
situation regarding initial probability distribu- 
tions in momentum and final probability distri- 
butions in space. In standard slit diffraction ex- 
periments, the momentum distribution of parti- 
cles just behind the slits is probed by the spa- 
tial distribution of the particles at the detectors 
very far behind from slits. Once this close anal- 
ogy between diffraction experiments and recent 
Bose condensed experiments is appreciated, it 
then becomes evident why these measurements 
have exhibited interference effects on the meso- 
scopic scale of the trapped Bose condensates. 

Our brief review of diffraction theory from slits 
is as follows: (i) A particle with high velocity v 
is normally incident on a screen with one or more 
slits. The slits have a length scale d and the high 
incoming momentum is defined by the inequality 
Mv >> (Ti/d). (ii) After a time period 



t s w (D/v) (slit to detector time), 
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the particle hits some detector at a distance D 
behind the slits, (iii) Before the particle hits the 
slits from incidence side of the screen, the mo- 
mentum parallel to the screen is zero. If the par- 
ticle is squeezed through the slits to the emission 
side of the screen at time zero, then the quan- 
tum mechanical uncertainty principle forces the 
momentum amplitude 



A( P ) 
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The momentum component p is in the re- 
direction, and lives in the plane of the screen with 
the slits. In that plane, p is normal to the long 
slit direction, (iv) Far behind the slits, the spatial 
distribution in x near the detectors is given by a 
one-dimensional version of Eq.(5); It is 



|#M.)| 2 « / \A{p)\H(x~{pt s /M)) 
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From Eqs.(10) and (12) follows the general slit 
diffraction formula 



For example, for diffraction through one slit hav- 
ing a width w, Eq.(ll) reads 



1 
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so that Eq.(13) yields the one slit diffraction pat- 
tern 
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For two slits, each of width w and separated by 
a distance b which obeys b > w, the momentum 
amplitudes obey 



A {2 - S ut)(p) = V2 cos (|jr)A 



(l-slit) 
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leading to the diffraction pattern 

(Ad vbx 
^-^)|^(i-^ t) (x)| 2 .(17) 

Thus, one sees that the momentum distribution 
near the slits give rise to the diffraction pattern in 
space near the detectors, via Eq.(13). Let us now 
turn to the quantum interference diffraction pat- 
tern expected from the momentum distribution 
of a trapped Boson condensate. 

In the zero temperature Gross pl[|- 
Pitaevskii|l5| theory for dilute Bosons in a trap 
potential U(r), the variational trial ground state 
is taken to be "quantum coherent" in a sense that 
is closely analogous to coherent states in quantum 
optics. The trial ground state is a normalized 
eigenstate of the Boson field operator, 



7/>(r)|>3> >= \&(r)|* > 



(18) 



The Bose condensation order parameter \&(r) is 
then chosen to minimize the energy functional 



£ =< ^\H-fj,N\^ > 



(19) 



where [i is the chemical potential. Eq.(19) is eval- 
uated using the two body scattering length a; 

8 = J <jV(/i - /i)* + (27r?i 2 a/M)|*| 4 }d 3 r, (20) 

where 
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Minimizing Eq.(20) yields 

/i*(r) + (47ra 2 a|*(r)| 2 /Af)*(r) = /z#(r). (22) 

In the trapped Boson experiments, one starts 
the design of the trap with an inhomogeneous 
magnetic field so that Uo(r) behaves as an 
anisotropic oscillator potential with a single min- 
imum at the origin. Such a single minimum in 
Uo(r) yields for the order parameter v I / o(r) a sin- 
gle maximum. This single maximum in \&o( r ) 
would behave as single source for the conden- 
sate explosion once the trap potential is removed. 
However, the single minimum in Uq(t) (i.e. the 
single maximum in 'l'o(r)) system is hard to 
build because the integer oriented spin atoms slip 
through a "hole" at the trapping potential mini- 
mum. 

To fix this potential, some workers have 
plugged up the hole at the bottom of the mag- 
netic trap with a LASER beam. The plug pro- 
duces (as a side effect) an effective U (r) with two 
minima. Thus, traps with a plug can be (and 
have been) designed || so that the order param- 
eter vE'(r) has two maxima. In such a case, when 
the trap is removed and the condensate explodes 
away from the trap region, there will be a dou- 
ble source for the condensate Bosons. If these 
two sources are coherent, then there was in the 
laboratory amplitude interference yielding order 
parameter fringes. These fringes were known to 
be analogous to the interference fringes observed 
in conventional "two slit diffraction" . If the two 
sources were incoherent, then such interference 
fringes would not have occured. 

It is to be expected that Bose condensed sys- 
tems (even at finite temperature) will exhibit a 
well defined coherent order parameter \P(r) whose 
condensate momentum distribution 



iVo(p) 
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will properly reflect the two maxima in ^(r). 
Such a double maxima describing two coherent 
sources duplicate in p(r, t) the interference pat- 
tern closely analogous to two slit diffraction. Af- 
ter the trapping potential is removed, the con- 
densate atom contribution to the density is given 



by 

p (r,t) « J N o (p)5(r-(pt/M))j0^ . (24) 

Reasoning as in Eqs.(16) and (17), and employing 
a separation vector b which describes the position 
of the second peak in ^(r) relative to the first 
peak in ^(r), we find that 

N ( 2 -peak) (P) ~ J(P,b)iVo(l —peak) 

where the interference function 

7(p,b) = (l+?7cos(p-b/7i)), (26) 

and < i] < 1 describes the "visibility" of 
the interference fringes when the two peaks (i.e. 
sources) are only partially coherent. 

Finally, we note that the oscillations in mo- 
mentum p implicit in Eq.(26) for the interfer- 
ence function i(p,b), have appeared as oscilla- 
tions in space r for a "one frame snap shot" || 
of po( r ,tfixed) fixed in time. Interference is also 
present as a time oscillation in po( r fixed, t) for a 
movie picture of atoms in the neighborhood of a 
fixed point in space Eg]. 

The point is that the momentum, position 
and time are related by the "classical path rule" 
r = vt = (pt/M), but with a weighted momen- 
tum distribution. This essence of the diffraction 
of particle beams as described by quantum me- 
chanical amplitude superposition, should also be 
possible to observe in systems for which the one 
rotates the potential about an axis displaced from 
the symmetry axis. Mesoscopic quantum interef- 
erence should also be present in TOP traps 0. 
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